Equivalent qubit dynamics under classical and quantum noise 
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We study the dynamics of quantum systems under classical and quantum noise, focusing on de- 
coherence in qubit systems. Classical noise is described by a random process leading to a stochastic 
temporal evolution of a closed quantum system, whereas quantum noise originates from the coupling 
of the microscopic quantum system to its macroscopic environment. We derive deterministic master 
equations describing the average evolution of the quantum system under classical continuous-time 
Markovian noise and two sets of master equations under quantum noise. Strikingly, these three 
equations of motion are shown to be equivalent in the case of classical random telegraph noise 
and proper quantum environments. Hence fully quantum-mechanical models within the Born ap- 
proximation can be mapped to a quantum system under classical noise. Furthermore, we apply 
the derived equations together with pulse optimization techniques to achieve high-fidelity one-qubit 
operations under random telegraph noise, and hence fight decoherence in these systems of great 
practical interest. 



I. INTRODUCTION 

A quantum computer is an emerging computational de- 
vice superior to its classical counterpart in certain tasks 
of great practical significance, e.g., factorization of large 
integers unsorted database search and the simu- 
lation of quantum systems [1, |3| ■ The main concern in 
building a working large-scale quantum computer is the 
undesired coupling of quantum bits, qubits, to their en- 
vironment. This coupling results in errors and loss of 
purity, a phenomenon generally referred to as decoher- 
ence. However, an important branch of quantum infor- 
mation science, fault tolerant quantum computing Q, 
states that if the qubit operations can be achieved with 
an error probability below a certain constant threshold, 
the total computational error can actually be rendered 
arbitrarily small. Thus studies on the effects of noise and 
decoherence on the qubit dynamics and especially direct 
suppression of the induced errors in qubit operations are 
of great interest. 

Recently, dynamical decoupling of quantum systems 
from their environments has been under active research. 
In the original scenario [6l|, hard, infinitely fast control 
pulses were applied to the qubit to effectively turn off 
the coupling to its environment, and hence to preserve 
its state by suppressing decoherence. It has been shown 
that it is also possible to generate qubit operations with- 
out disturbing the decoupling process [7] and that fast 
soft pulses suffice for the control [8|]. Furthermore, the 
effectiveness of dynamical decoupling has been improved 
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by introducing randomness to the control [9|, llO| . On the 
other hand, a direct pulse optimization method to obtain 
high-fidelity qubit operations in the presence of classical 
random telegraph noise (RTN) in the qubit energy split- 
ting was reported in Ref. 11]. 

In this paper, we derive deterministic master equations 
for the average temporal evolution of a quantum sys- 
tem under classical continuous-time Markovian noise. In 
the case of a qubit under RTN, our equations reduce to 
two coupled one-qubit master equations. This formula- 
tion allows for great speed-up in dynamical decoupling as 
well as direct pulse optimization schemes for high-fidelity 
quantum gates compared with ave rag ing over many noise 
realizations as was done in Ref. [ll|. Since the master 
equations are not stochastic, they are also well suited for 
algebraic manipulations. 

For quantum noise, we study the dynamics of a sys- 
tem coupled to a macroscopic quantum environment 
within the standard Born approximation. In addition, 
we present another set of master equations for a sys- 
tem under noise due to coupling to an impurity which is 
furthermore coupled to a Markovian environment. The 
latter case has been extensively studied, for example, in 
Refs. for suppression of random telegraph and 

1// noise in solid state qubits. Even though the obtained 
master equations for quantum noise are already known, 
we demonstrate an interesting connection between the 
different approaches — the master equations for classical 
noise are equivalent to those obtained in the case of quan- 
tum noise for certain couplings of the system to the envi- 
ronment. This result shows how classical noise can arise 
from coupling to a quantum bath, and that in some cases 
classical noise can accurately model quantum noise. We 
stress that our master equations yield the complete dy- 
namics of the quantum system under noise in the Born 
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approximation, not only for example damping rates stud- 
ied in Ref. [H. 



III. OPEN QUANTUM SYSTEM IN THE BORN 
APPROXIMATION 



II. CLOSED QUANTUM SYSTEM UNDER 
CLASSICAL MARKOVIAN NOISE 



We begin by considering a quantum system influenced 
by classical continuous-time Markovian noise with N dis- 
crete states, each corresponding to a different system 
Hamiltonian. The average outcome of any measurement 
applied to the system is obtained from p{t), the density 
operator averaged over all possible noise trajectories. Al- 
though each noise trajectory alone corresponds to unitary 
evolution, the averaged evolution can be nonunitary. 

The Markovian assumption implies that the noise has 
no memory, i.e., at time t the state of the noise process 
is fully described by the probabilities Pk(t) of the noise 
states, and the probability flow is governed by the tran- 
sition rate matrix 7 — {'^kj} as 



N 



(1) 



where X)fc^fc(^) ~ 1- We define the operator pk{t) as 
the density operator averaged over all noise trajectories 
occupying the state k at time instant t, and normalized 
such that Tr[yOfc(t)] = Pk{t). The total state operator of 
the system is obtained as 



N 



(2) 



The evolution of each pk{t) under an infinitesimal time 
interval is given by the sum of all possible quantum evo- 
lutions weighted by their probabilities. Thus we obtain 
a closed set of ordinary master equations for the condi- 
tional density operators Pfc(i): 



N 



dtPk{t) 



[Hk{t),pk{t)]+Y.^k,pAt)^ (3) 



where Hk (t) is the system Hamiltonian corresponding to 
the fcth noise state. One can verify the consistency of 
Eqs. H]) and ^ by taking a trace of Eq. ([3]), which re- 
sults in Eq. JT]) . We point out that the conditional density 
operators pk (t) are introduced only for calculational pur- 
poses to obtain the average density operator p{t) accord- 
ing to Eq. For RTN, i.e., a single bistable fluctuator, 
Eqs. ([3]) reduce to 

dtP±{t) [H±it),p±it)] ± - (4) 

where the correlation time of the noise is denoted by Tc. 
Classical RTN similar to Eq. ([4]) has been previously 
studied for example in the case of macroscopic quantum 
tunneling in Josephson devices [l^ . We apply Eq. ^ to 
finite-dimensional qubit systems in Sec. |Vl 



Let us turn our attention to quantum noise by study- 
ing a quantum system coupled to its environment. The 
temporal evolution of the total density operator pt is gov- 
erned by the Hamiltonian 

H = Hs(S)Ic+Is(S He + gK^ g, (5) 

where the pure system and environmental Hamiltonians 
are denoted by Hg and 77c, respectively, and the inter- 
action of the system with the environment is described 
by Hint ■ Without loss of generality, we can assume that 
the expectation value of the operator E over the envi- 
ronmental degrees of freedom vanishes. To obtain the 
lowest-order corrections in the coupling strength g to the 
temporal evolution of the system, we write the master 
equation for pt in the interaction picture with respect to 
Ho as dtpt = ^[Hint,Pt], where pt := Ul{t)ptUo{t) and 
Hint ■= UQ{t)HintUo{t). The evolution operator Ua{t) is 
expressed using the time ordering operator T as Uoit) — 
Tew /o ■'^o(^)'i*. The master equation can be expressed 
equivalently in an integro-differential form as [l^ 



dtPt = 4 Hintit),pt{0) 
iTi 



(6) 



{ihf 



Hint{t), Hint{s),Pt{s) 



ds. 



We proceed by assuming that the coupling of the sys- 
tem to the environment is sufficiently weak for the Born 
approximation to be valid, i.e., pt ~ ps ® Pc- Thus the 
first term on the right hand side of Eq. ([6|), containing 
operator E only to the first order, vanishes in taking the 
trace over the environmental degrees of freedom. The 
remaining part of the equation yields 



dtPs 



g 



[[K(t),K{s)p,{s] 



k{t),-p,{s)K{s) C{Q,t 



C(t-s,0) 
s)}ds, (7) 



where the forward and backward autocorrelation func- 
tions are defined as C{t - s, 0) = Tr{£^(t - s)E{{))pc} and 
C(0, t - s) = TT{E{{))E{t - s)pc}- 

Since the operators E at different instants of time do 
not necessarily commute, Eq. ([7]) cannot in general be 
simplified. However, for environments at high enough 
temperatures T the Kubo-Martin-Schwinger boundary 
condition C{t - s,0) = C[0,t - s + ih/{kBT)] yields 
C(i - s, 0) « C(0, s) =: Cc{t - s), and hence Eq. ^ 
reduces to 



dtPs = 



g 



c,{t-s) k{t), k{s),ps{s) 



(8) 



For example, quantum noise arising from a trapping cen- 
ter hybridized with the Fermi sea in a superconductor 
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at temperatures much higher than the Unewidth of the 
trap results in RTN |18| with the autocorrelation function 

2 |a-t| 

Cc{t — s)—e . Surprisingly, insertion of this corre- 

lation function into Eq. ([8]) results in the same temporal 
evolution of the average density operator as the master 
equations (|3|) for classical RTN. This equality can be for- 
mally verified by showing that if Eq. ([8]) holds, 



P± 



2^' 2ih 



K{s),Ms) ds (9) 



is a solution to Eq. ([4]) with H± = ±gK, and that 
p+ + P- ~ Ps- Thus we have shown that quantum noise 
arising from a realistic environment can be accurately 
modeled with classical noise within the Born approxima- 
tion. We note that although Eq. ([8]) seemingly includes 
memory effects, i.e., the time derivative of the density 
matrix at time instant t depends on the density matrix 
at earlier times, the above proved equivalence shows that 
no memory is required in the case of RTN if the dimen- 
sion of the differential equation system is doubled. 

The above approach can also be employed when the 
influence of the environment is modeled classically by 
including a stochastic noise term in the system Hamil- 
tonian. A stochastic system Hamiltonian analogous to 
the Hamiltonian in Eq. ([5]) reads H = Hs + g£^{t)K, 
where £,{t) is determined from some stochastic process. 
The resulting equations are similar to those presented 
above, except that each expectation value over the envi- 
ronmental degrees of freedom is to be replaced with an 
ensemble average over noise trajectories. In particular, 
Eq. ([8]) retains its form. 



IV. OPEN QUANTUM SYSTEM WITH 
LINDBLAD DAMPING 



Above, we have shown that a classical noise source 
can be used as an effective description of a decoher- 
ence process arising from the microscopic dynamics of 
a many-body quantum system within the Born approx- 
imation. Below, we consider a physically relevant deco- 
herence source of a superconducting qubit — a defect cou- 
pled to a Markovian environment [isl . [l^ , and show 
by tracing out the environmental degrees of freedom how 
this model is also mapped analytically to Eq. (g]) de- 
scribing the average system dynamics under RTN. This 
central result motivates and further justifies the use of 
classical noise in master equations of quantum dynamics. 
The reduction of the complicated many-body dynamics 
to the form presented above allows significant conceptual 
and practical simplification of the problem. 

Decoherence mediated by fermionic vacancies in a 
Markovian environment has been studied in Refs. p^.[l3l. 
[T3 | as a possible source of RTN and 1// noise in Joseph- 
son qubits. In addition, an equivalent defect has been 
shown to result in the RTN spectrum [l^. Motivated 
by these results, we consider a microscopic quantum sys- 
tem coupled to a defect through the projection operators 



on the eigenbasis of the defect. Without loss of general- 
ity, the total Hamiltonian of this system and an electron 
band coupled to the defect can be expressed as 



db, 



(10) 



where the operator Kg of the microscopic quantum sys- 
tem is responsible for the coupling to the defect, Uzd is the 
Pauli z matrix operating in the eigenbasis of the defect, e 
is the energy splitting of the defect, and Hdh couples the 
defect to an electron band represented by the Hamilto- 
nian i/b- The defect is essentially a bistable quantum 
fiuctuator which can also be represented by a fermionic 
vacancy with a population varying between and 1 as 
electrons tunnel into and out of the band. By eliminat- 
ing the band electrons from the dynamics of the defect 
in the Born-Markov approximation, we obtain effective 
terms operating on the defect part of the density operator 
space represented by the Lindblad operator 

C-d{pd} = (2o-"pdCr+ - cr+O-^pd - pdCr+fT") 

+ — (2CT+PdO'" - (J^CF^Pd - Pd<y^<7^) ■ (11) 

The first three terms describe the process of electron 
tunneling out of the defect with rate Fi and the last 
three terms describe the inverse process with rate 
The rates obey the detailed balance condition ri/r2 — 
exp[e/ (fceT)], where T is the temperature of the electron 
band. The equation of motion for the density opera- 
tor pds of the microscopic system and the defect can be 
written as 



dtPds =-T[Id ®Hs + ^(T,_d 

in 2 
+ (£d ®/s){pdJ. 



' h + CTzd ® Ks, Pds 



By expressing the density operator pds in the form 



Pds 



P++ P+- 
P-+ P-- 



(12) 



(13) 



and inserting it into Eq. ()12|) . we obtain 



dtPds =- 



l_([H,,p++] [H,,p+4\ e 
th \[Hs,p-+] [Hs,p—]J ^ 2ih 

If [K,,p++] {K,,p+-} 
^ ih \~{K,,P-+} -[K,,P—] 

El /'-2p++ -p+-\,Elf^P- 
^ 2 I -P-+ 2p++ ; ^ 2 [-P- 



2p+_ 
2p-+ 



(14) 



-P+- 
-2p— 



Equation p4|) shows that the diagonal blocks of pds de- 
couple from the off-diagonal ones. The diagonal blocks 
are of primary importance since the density operator 
of the system can be expressed as ps = Trd{pds} — 
p^^ + p , and hence its dynamics is completely de- 
termined by the two coupled master equations with the 
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same dimension as ps- In fact, the dynamics of pds in the 
high temperature hmit corresponds to Eq. which is 
observed by denoting p±± := p±, Ti — T2 1/tc, and 
H± := Hs ± Ks- Thus we have shown that the decoher- 
ence arising from the couphng of the system to the defect 
can be modeled with classical RTN. For finite tempera- 
tures, Eq. can be recast into the form of Eqs. ^ 
showing that the equivalence of classical and quantum 
dynamics does not necessarily arise due to the infinite 
temperature limit. 



V. HIGH-FIDELITY ONE-QUBIT NOT GATES 

To demonstrate the computational effectiveness of the 
derived master equations, we implement high-fid elity 
quantum gates for a qubit under RTN as in Ref. 
We assume that the qubit dynamics can be controlled in 
the (Ta, direction and that the noise acts in the az direc- 
tion with the strength A, i.e., H± — [a{t)ax ± Aaz]/2 in 
Eq. The strength of the control field a{t) is assumed 
to be bounded by Cmax- 



(a) 




10 20 30 40 

Tc/(?i/amax) 



50 




10 20 30 40 

Tc/ (?i/amax) 



50 



FIG. 1: Fidelities for the NOT gate as functions of the correla- 
tion time Tc for a 7r-pulse (dotted line), CORPSE (dash-dotted 
line), short CORPSE (dashed line), and gradient optimization 
(solid line). The RTN strength A is set to (a) 0.125 x flmax 
and (b) 0.25 x flmax- The gate fidelity for a gate U is defined as 
{Tt{U poW p}) , where the average is calculated over all pure 
initial states po, see Ref. 

Figure 1 shows gate fidelities $ for the NOT gate, i.e.. 
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FIG. 2: Optimized pulse sequences yielding the highest gate 
fidelities for correlation times (a) 5h/ama,x, (b) 20?J,/amax, and 
(c) SOR/omax- The strength of the noise is chosen to be A = 
0.125aniax corresponding to Fig. [T]^a). 



the unitary operator cTa,, as functions of the noise cor- 
relation time, obtained using composite pulse sequences: 
TT-pulse, compensation of off-resonance with a pulse se- 
quence (CORPSE), and short CORPSE To 
achieve optimized gate fidelity, we employed a variant 
of the gradient ascent pulse engineering method [l^] of 
a piecewise constant function a(t). Figure 2 shows op- 
timized pulse sequences for different correlation times of 
the RTN fluctuator. We observe that for short corre- 
lation times the optimized sequence is very close to a 
TT-pulse and for long correlation times the pulse sequence 
assumes a shape similar to short CORPSE. 

In contrast to similar results presented in Ref. for 
A = 0.125 X Umax, the curves in Fig. [1] do not show 
the statistical errors arising from the finite sampling of 
the noise. One of the conclusions in Ref. fTT] was that 
the gradient optimization yields only a marginal improve- 
ment to the gate fidelity over the most efficient composite 
pulse sequence, see Fig. [TJa). However, we find that this 
conclusion is valid only for very weak noise, and for ex- 
ample the noise strength A = 0.25 x Omax employed in 
Fig. [Hb) suffices to render gradient optimization clearly 
the most efficient method considered to fight the deco- 
herence in this system. 



VI. CONCLUSION 

We have shown how classical random telegraph noise 
can arise directly from coupling of a quantum system to a 
macroscopic quantum environment or indirectly through 
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an impurity coupled to a Markovian environment. The 
average dynamics of the quantum systems were found to 
coincide under classical and quantum noise within the 
Born approximation. Our observations justify the uti- 
lization of classical noise as a model for quantum noise 
for certain systems. Furthermore, the results presented 
in this paper introduce a possibility of high-performance 
pulse optimization for qubits in noisy environments, and 
offer an interesting point of view for studies on the effect 
of noise on quantum systems in general. In the future, 
we concentrate on the optimization of multi-qubit opera- 
tions under noise, study the effects of different quantum 
baths on the qubit dynamics, and aim to generalize the 
obtained equations for continuous and more complicated 
noise models. The obtained results for the average quan- 
tum dynamics under noise are not restricted to qubits 



or their multi-level equivalents, qudits, and hence can be 
employed in studying, e.g., dilute Bose- Einstein conden- 
sates in the noninteracting limit. 
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